Ground state properties of one-dimensional ultracold Bose gases in a hard-wall trap 
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We investigate the ground state of the system of N bosons enclosed in a hard-wall trap interacting 
via a repulsive or attractive 5-function potential. Based on the Bethe ansatz method, the explicit 
ground state wave function is derived and the corresponding Bethe ansatz equations are solved 
numerically for the full physical regime from the Tonks limit to the strongly attractive limit. It 
is shown that the solution takes different form in different regime. We also evaluate the one body 
density matrix and second-order correlation function of the ground state for finite systems. In the 
Tonks limit the density profiles display the Fermi-like behavior, while in the strongly attractive limit 
the Bosons form a bound state of N atoms corresponding to the iV-string solution. The density 
profiles show the continuous crossover behavior in the entire regime. Further the correlation function 
indicates that the Bose atoms bunch closer as the interaction constant decreases. 

PACS numbers: 03.75.Hh,05.30. Jp,67.40.Db,03.65.-w 



The physics of one-dimensional (ID) cold atoms has 
recently attracted a great amount of attention due to 
tremendous experimental progress in the realization of 
trapped ID cold atom systems 0, H 0. A ID 
quantum gas is obtained by tightly confining the particle 
motion in two directions to zero point oscillations [|| . As 
the radial degrees of freedom is frozen, the quantum gas 
is effectively described by a ID model along the longitu- 
dinal direction 0. Experimentally, a ID Bose gas can 
be realized either by means of anisotropic magnetic trap 
or two-dimensional optical lattice potentials. In paral- 
lel, the wide exploitation of the Feshbach resonance to 
control the scattering length of the atoms allowed exper- 
imental access to the full regime of interactions both from 
weakly interacting limit to strongly interacting limit and 
from the repulsive regime to the attractive regime by sim- 
ply tuning a magnetic field. Very recently, several groups 
have reported the observation of a ID Tonks-Girardeau 
(TG) gas 0, 0, 13 j which provides a textbook example 
where atom-atom interactions play a critical role and the 
mean- field theory fails to obtain reasonable results 0. 
Theoretically, the effect of dimensionality and correla- 
tion effect of bosonic system have been investigated ex- 
tensively [l|![l3[ll|HHQHE[3[l3andis 
being paid more and more attention. 

It is well known that there exist some exactly solved 
ID interacting models HJ El |H E3, however, 
are not directly applicable to the system trapped in 
an external harmonic potential and the Gross-Pitaviskii 
(GP) theory is widely adopted in dealing with the sys- 
tem of Bose-Einstein condensations (BECs) with weak 
interactions. In spite of its great success in account- 
ing for the basic experimental observations, the GP 
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equation is essentially based on a mean-field approx- 
imation and suffers from various shortcomings, espe- 
cially when applied to the strongly correlated systems. 
Some recent works have shown the limitations of GP 
theory in the density distribution of ID trapped gas 
in strongly interacting regime jlfij. the overestimation of 
interference 0|, and the instability of GP equation with 
attractive interactions [24|. For a trapped system, the 
attractive interactions may compensate the kinetic en- 
ergy of the condensate and lead to a stable soliton solu- 
tion. Recently, the ID GP equation under box boundary 
condition was solved analytically for both repulsive and 
attractive condensates [25j . On the other hand, the ex- 
perimental progress in building up square well trap[2rjj 
and optical box trap j^, gives rise to the hope to di- 
rectly study the physics in the textbook geometry of a 
"particle in a box" . It is not clear whether the solution of 
GP equation is good enough to describe the interacting 
Bose gas in the hard-wall trap, especially for the attrac- 
tive case and under the strongly interacting limit. Fortu- 
nately, in a hard-wall trap, the corresponding interacting 
model is integrable and its exact solution has been ob- 
tained with the Bethe ansatz method for the repulsive 
interaction in a seminal paper by Gaudinpsj. So far, 
there has been a growing interest in the exactly solved 
models in the hard- wall trap [2^, |3(], E3- In spite of 
the long history of the integrable Bose interacting model 
(Lieb-Liniger model), the case with attractive interaction 
draws less attention|23, H3, H3 and most of the studies 
have focused on the ground state energy for the periodic 
boundary systems. 

In this paper, we investigate the density distribution of 
the ground state of the ID Bose gases in an infinite deep 
square potential well. Different from the system with pe- 
riodic boundary condition where the density distribution 
is a constant and is irrelevant to the strength of interac- 
tion, the density profile for a trapped bosonic system is 
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FIG. 1: (color online) The density of state in quasi- 
momentum k space for the ground state for N = 200 and 
c = 0.1 (solid line), c = 1 (dashed lines), c = 10 (dash dot 
lines), c = 100 (dash dot dot lines). Inset: The numerical 
solutions of the Bethe ansatz equations for N = 1000 and 
c= 10. 



shown to be sensitive to its interaction. The dependence 
of one body density matrix on the interaction between 
atoms and the effects of finite size are studied for both 
the repulsive and attractive interactions. The wave func- 
tion of ground state is obtained by numerically solving 
the set of Bethe ansatz equations. 

We consider N particles with the S interaction in one 
dimensional box with length L. The Schrodinger equa- 
tion can be formulated as (the natural unit is used) 



N 

£ 

i=l 



dxj 



2c Y 6 ^ 

l<i<j<N 



where 2c is the interaction strength between atoms re- 
lated with the s -wave scattering length 0, E3 , which 
can be tuned from — oo to +oo by the Feshbach resonance 
or confinement induced resonance. When the interaction 
between atoms is repulsive, c > 0, and in the attractive 
case c < 0. The important parameter characterizing dif- 
ferent physical regimes of ID quantum gas is 7 = cp , 
where p = N/L. We will study the full physical regimes 
—00 < c < 00. The wavefunction can be written as 

*(xi,--- ,x N ) = ^2<p(x pi> x pa ,--- ,x PN ) 



x 9 (x Pl < x P2 ■ ■ ■ < x Pn ) , 



in which 
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ViiVii ' ' ' iPn is one of the permutations of 1, • ■ • , JV, and 
is the sum of all permutations. 6{x — y) is the step 



function. The wavefunction tp (x Pl , x p 



1 x p 



could 



be obtained by the permutation of ip(x±,--- ,xn) ac- 
cording to the symmetry condition of the Boson wave 
function. It turns out that the original problem is equiv- 
alent to solving the equation 
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2c 



l<i<j<N 



S (Xi - Xj) 



,X N ) 



tp(xi,--- ,X N ) 

(1) 



in the region < X\ < X2 < ■ ■ ■ < xn < L. We take the 
wave function tp (x±, ■ ■ ■ ,xn) as the Bethe ansatz type 



Ap exp i Tjk p 



ip(x 1) x 2 , ■ ■ ■ ,x N ) = Y 

P,r 1: ...,r N 

(2) 

where rj = ± indicates that the particles move toward 
the right or the left. Substituting eq.J2J) into eq.QJ and 
using the open boundary conditions 



<p (0,2:2, • • • , Xn) = <P (X1,X 2 

we have the Bethe ansatz equations 



,L) = 0, 



exp (i2kjL) 



n 



iki + ikj — c iki — ikj + c 
iki + ikj + c iki — ikj — c 



(3) 



with j = 1,2, ■ ■ ■ , N. The energy eigenvalue is E — 

SjLi hj an( i the total momentum is k = YljLi kj- 

Taking the logarithm of Bethe ansatz equations, we 
have 



arctan ■ 



N f c 

kjL = rij-K + > I arctan — 

^ \ kj + ki 

(4) 

with j = 1, • • • ,N. Here {rij} is a set of integer which 
determines an eigenstate and for the ground state rij = 1 
(1 < j < N). Alternatively, eq. (J1J has the form of 



kjL = rijir 
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E 



/ kj + k[ 
arctan — 

V c 



arctan ■ 



... (5) 
The above two equations are consistent but the choice of 

{n'j\ is different from that of {rij}. For the latter case, 

we should have n'j = j(l<j<N) for the ground state. 

The {kj}, and thus the wave function can be decided by 

the set of transcendental equations e q s. fit - In general, 

the eqs.© are used in the literature |2lll22l | for the model 

with repulsive interaction, because it is more convenient 

to extend it to deal with the problem of thermodynamics. 

For c > 0, the set of fcj is unique and real, therefore the 

density of state in k space p(k) can be formulated as |2(j 



Lp 



kj + kj + i 



(6) 
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However, such a definition does not make sense when 
c < 0. In the later, through an example with N = 3, 
we will show that in the attractive case the solution is 
not unique corresponding to a given set of {rij} and the 
ground state is decided by comparing the energy eigen- 
values of different solutions. 

With some algebraic calculation, the wave function has 
the following explicit form 



tp(xi,x 2 , ■ ■ ■ ,x N ) 
Ap £p exp 



p 



JV-l 
Kj 



exp (ik PN L) sin (k Pl xi) 



Kj<N \ Kj 
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FIG. 2: (color online) The density profiles p(x) of the ground 
state for different repulsive interaction constants and N = 4. 



w Q h = arctan ■ 



arctan ■ 



and 



N 

A PiP2 ... Pn = Y[ (ikp, — ik Pj + c) (ik Pl + ik Pj + c) . 

Here G p denotes a +(— ) sign factor associated with even 
(odd) permutations. In terms of the ground state wave 
function '5 (x\, • • ■ , Xjy), the one body density matrix is 
defined as 



p(x) 



N J £ dx 2 ■ --dx N |\fr (x,x 2 ,x 3 , ■ ■ ■ ,x N )\* 
J dxi ■ --dx N |* (xi,x 2 , ■ ■ ■ ,x N )\ 2 



In the following calculation, we will let L = 1 through 
the paper. Firstly, the repulsive interaction case is con- 
sidered. The density of ground state in quasi-momentum 
k space for different interaction constants is plotted in 
Fig. 1. It is shown that the density is suppressed for k 
close to zero because of the confinement by the infinite- 
depth well. The numerical results of eqs. 10} are given 
for N = 1000 and c = 10 in the inset. Bethe ansatz 
equations uniquely decide the value of fcj, the wave func- 
tion of the ground state and the one body density matrix 
of the system. In Fig. 2, we show the one body den- 
sity matrix of the ground state for N = 4 for different 
interaction constant. When there is no interaction be- 
tween atoms (c = 0), all of the atoms have the same 
quasi momentum, kj — it (1 < j < N), corresponding to 
condensation of the ideal Bose gas. The one body den- 
sity matrix is the sum of the density of N independent 
bosonic atoms lying in the ground state. As the interac- 
tion increases further, the half-width of the density of the 
system becomes larger and larger with increasing inter- 
action. At c — 10, the density distribution shows already 
Fermion-like behavior. As the further increase of interac- 
tion the Bosons display the same density profile as that of 




FIG. 3: Three sets of solution (squares, triangles, circles) for 
3 Bosons in a hard-wall trap with the attractive interaction 
c = — 1. As comparison, diamonds correspond to the ground 
state solution of a periodic 3-Boson system on a ring of length 
L = 1 for c = -1. 



N noninteracting spinless Fermions and we have kj ~ jtt 
(1 < j < N) at c ~ 1000. 

In the case of attractive interaction, c < 0, eqs. (@J 
have solutions with complex quasi-momentum. The case 
for periodical boundary condition has been investigated 
in reference [2^, |3^, |3J| . Due to the zero-point energy 
for a confined system, the solution corresponding to the 
ground state is not purely imaginary like that for the 
Bethe ansatz equations with periodical boundary condi- 
tion. For the two body problem, the complex solutions 
of the Bethe ansatz equations take the form of the two- 
string solution 

k\ = a — iA, 
k 2 = a + iA, 

where a and A are real and 2a is the momentum of the 
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system. In general, by an iV-string solution we mean a 
solution where N momenta k{ possess the same real part. 
Then the eqs. (@J have the following form 



n 1 +n 2 c 

aL = 7r + arctan — , 

2 2a 

1 , 2A - c 
AL = -In— , 

2 2A + c 



(7) 



where ni = n% = 1 for the ground state. The correspond- 
ing energy of the system is E = 2a 2 — 2 A 2 . In terms of 
the momentum k = 2a, the energy can be represented 
as E = k 2 /2 — 2A 2 , which implies that the state corre- 
sponding to a two-string solution can be regarded as a 
bound state of two atoms with a binding energy —2 A 2 
and a doubled mass. For convenience, we also refer to 
the bound state corresponding to a two-string solution 
as a dimer state in the following text. 

For the three-body problem, we assume that the solu- 
tion has the following form: k% — a — iA, ki — a + iA 
and &3 = (3, therefore the Bethe ansatz equations can 
be represented with a, (3 and A. For the case with pe- 
riodic boundary condition, there exists a 3-string solu- 
tion corresponding to a = /3, i.e. the trimer state (the 
bound state of three particles) |22j, |33j ■ But we can prove 
that no 3-string solution could be formed in a box for fi- 
nite attractive interaction (see the appendix for details). 
That means the open boundary condition or a confine- 
ment tends to prevent the formation of a trimer state. 
We show the solutions of three-body problem for period- 
ical boundary condition (Diamonds) and open boundary 
condition (Squares, Triangles and Circles denoting three 
sets of different solutions) for c = — 1 in Fig. 3. It should 
be noticed that the set of solution denoted with squares 
is not an exact 3-string solution although a ~ /3, which is 
a dimer plus a single atom. By comparing E, the trimer- 
like solution (squares) has the lowest energy eigenvalue 
and thus it represents the ground state of the system. 

For the system with four atoms, the general complex 
solutions can be assumed as k\ = a — iAi, ki = a + iAi, 
&3 = j3 — iA2 and Lt = (3 + iA2- By solving the 
Bethe ansatz equations numerically, it turns out that the 
ground state corresponds to a string solution of length 
four (a = (3 but Ai 7^ A 2 ) for the large attractive interac- 
tion, whereas the ground state is a scattering state of two 
dimers for the weak attractive interaction. For N > 4 the 
numerical results show that generally the ground state 
of the system in the weak attractive regime is a scat- 
tering state of N/2 dimers (N is even) or (N — l)/2 
dimers plus a single atom (N is odd). The state of N/2 
dimers is given in terms of N/2 two-string solutions. As 
the attractive interaction increases further, the ground 
state smoothly evolves into an intermediate regime char- 
acterized by an (N — M)-string solution plus M/2 two- 
string solutions with 2 < M < N, and finally it falls 
into the strong coupling regime corresponding to the N- 
string solution. To give an explicit example, we display 
the ground state solutions to the eqs.QJ for N = 10 
with different finite attractive interactions in Fig. 4. 
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FIG. 4: The solutions to the Bethe ansatz equations for N = 
10 with different attractive interaction constants, c takes the 
value of -10, -4.0, -3.0, -2.0, -1.8, -1.6, -1.4, -1.2, -1.1, 
-1.0, -0.9, -0.8, -0.7, -0.6, -0.5 (from left to right). Inset: 
(a) c = -0.8; (b) c takes the values of -0.4, -0.3, -0.2, -0.1 
(from left to right). 



Based on the numerical solutions of the Bethe ansatz 
equations, we find that 5 dimers form the ground state 
of the system when the attractive interaction is weak 
(—1.42 < c < 0). While in the limit of strong at- 
tractive interaction(c < —3.02), the ground state solu- 
tion is a 10-string solution. In the intermediate regime 
(—3.02 < c < —1.42), the ground state solution is com- 
posed of a 4-string solution plus three 2-string solutions 
for —1.72 < c < —1.42, a 6-string solution plus two 2- 
string solutions for —1.88 < c < —1.72, and an 8-string 
solution plus a 2-string solution for —3.02 < c < —1.88. 

For the periodic system, the A^-string ansatz is repre- 
sented as 



k j =a + i(N+l-2j)(c/2), j = l,- 



,N. 



The N-string ansatz is generally believed to be the solu- 
tion of the Bethe ansatz equations for c < in the limit 
L — > oo 22J. For the Bose gases in a finite-size hard- wall 
trap, we assume that the A^-string solution takes the fol- 
lowing form 



kj = a + i(N + 1 — 2j)(c/2 + Sj), 



3 = 1, 



,N, (8) 



where 8j is a set of small numbers and Sj — > as c — > 
— oo. Our numerical solutions to eqs.Q indicate that 
the solutions are precisely fitted by the A^-string ansatz 
of eqs.© in the limit of large attractive interaction. For 
instance, for the system with ten atoms the solution takes 
the values of kj = 0.498075-i(ll-2j)(5+£,-) with |^-| < 
10- 8 forc= -3.02 and/cj = 0.354213 - - 2j)(5 + Sj) 



with \5A < IQ- 13 for c 



TO. In the strong attractive 



limit, we can determine a analytically. With eq. JSJ and 
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the original Bethe ansatz equations the total momentum 
can be formulated as 



exp (i2kL) 



n 

l<j<N 



2a + i(jc) 



2a - i(jc) 



where k — Na and we have taken Si = 0. Taking the 
logarithm to the above equation, we have 



kL = 



ETC 
2 arctan — , 
2a 

l<j<N 



(9) 



It 



is clear that the 
oo and therefore 



where n = TV for the ground state 
momentum k is 7r in the limit of c 
a = tt/N in the strong attractive limit. 

In Fig. 5 we plot the one body density matrix of the 
ground state for TV = 2 (a) and TV = 4 (b). It is shown 
that as the attractive interaction increases, the central 
density of the system becomes large firstly and then less. 

In the limit of c -> -oo (c 1000 for TV = 2), the 

density profile matches the case of c = 0. The matching 
between them can be explained as a compounded particle 
with mass TVm located in its ground state, whose density 
distribution has the form of sin 2 x in the strong attractive 
limit. When the solution of Bethe ansatz equations take 
the form of eq.|[S}, the ground state energy can be rep- 
resented as E — k 2 /N - TV(TV 2 — l)c 2 /12, which implies 
that the state corresponding to the TV-string solution can 
be regarded as a bound state of TV atoms with a binding 
energy oc — c 2 and a mass TVm. It is also interesting to 
study the second order correlation function 



9i{x, c) 



N(N - 1) Jq dx 3 ■ ■ -dx N | VP (x,x,x 3 , ■ ■ ■ ,x N )\ 
Jq dxi ■ --dx N |W (xi,x 2 , ■ ■ ■ ,x N )\ 2 



In Fig. 6 we show the dependence of second-order cor- 
relation function upon the interaction. It indicates that 
the atoms tend to cluster together more easily for the 
attractive interaction and the atoms bunch closer as the 
interaction becomes stronger. For the repulsive interac- 
tions, the atoms avoid each other and the atom-bunching 
reduces and vanishes finally for increasing interactions, 
which is similar to the case for the periodical boundary 
condition (35| . 

In conclusion, by numerically solving the Bethe ansatz 
equations we investigate the ground state properties and 
obtain the density distribution function and the second- 
order correlation function of the ID Bose gases in a box 
of finite length L in all the physical regimes (— oo < c < 
+oo). In the limit of c — ► oo the Bose gas shows simi- 
lar behavior as that of noninteracting spinless Fermions, 
while in the limit of c — » — oo the Bose gas behaves as 
a compounded particle with mass Nm. In the case of 
weakly attractive interaction the ground state is com- 
posed of N/2 dimers (TV is even) or (TV — l)/2 dimers 
plus a single atom (TV is odd) . The second-order correla- 
tion function indicates that the atoms bunch closer as the 
interacting constant c decreases. Our results can cover 
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FIG. 5: (color online) The density profiles p(x) of the ground 
state for different attractive interaction constants, (a) TV = 2 
and (b) TV = 4. 



— 1- 


- 1 — I- 




— 1 


— 1 — 


-1 1 1 1 r 


1 • 






1.0 


1 




























t\T 


0.8 






























■J? 


0.6 














tT 


0.4 














■ 
















!!j 


0.2 














■ ra 


0.0 





















20 


40 60 8 
C 


100 




\ 



























FIG. 6: The second-order correlation function versus c for 
TV = 4. Inset: The second-order correlation function for re- 
pulsive interaction with TV = 4. 



the whole parameter regime beyond the mean field theory 
and display the continuous crossover behavior from the 
Tonks limit to the strongly attractive limit. Especially, 
the crossover behavior of the ground state for the case of 
attractive interactions is discussed in detail. Hopefully, 
our results based on the exact solution can provide a clear 
picture of the wavefunction for the BECs with attractive 
interaction in a trap and help us to gain some intuitive 
insight on the collapse of BEC. 
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APPENDIX A: THE STRING SOLUTION FOR 
THE ATTRACTIVE THREE-ATOM SYSTEM 

In this appendix, we study the three-particle system 
in detail. We firstly assume the 3-string solution for 3- 
atom system in the form of k\ = a — iA, ki = a + iA 
and &3 = a. In terms of a and A, then the Bethe ansatz 
equations take the following forms 



exp [(i2a 
i2a + A - 



2A) L] 

c i2a — c A 



c2A 



i2a + A + c i2a + c A + c2A + c 
exp [(i2a - 2A) L] 
i2a - c i2a - A- c2A + cA + c 
i2a + c i2a - A + c2A-cA-c' 
exp [i2aL] 

i2a + A - c i2a - A - c 



i2a + A + ci2a- A + c 



(Al) 



(A2) 



(A3) 



where we have used 4+£ 

A — c A+c 



1 in the right of the third 
equation. Substituting eq. (A3) into (Al) and (A2), we 
then get 



exp [2AL] 



i2a — A — c 



i2a 



■A 



■2A- 



i2a -A + c 



exp [-2AL] - 



i2a + A - c 



i2a + cA + c2A + c 
i2a - c2A + cA + c 



i2a + A + c i2a + c2A-cA-c 
The above two equations can be rewritten as 



exp [2AL] 



x exp 



/4a 2 + (A + c)' 
4a 2 + (A - cf 
A + c 



arctan ■ 



= exp 



2i arctan — 
2a 



arctan - 
2a 

A - c 2A - c 
A + c2A + c' 



A 



2a 



exp [-2AL] 



x exp 
exp 



Ua 2 + {A-cY 
Aa 2 + (A + cf 
A + c 



i arctan 

V 2a 



arctan ■ 



A 



2a 



2i arctan — 
2a 



2 A + c A + c 
2A-cA-c' 



Comparing the module or imaginary part of the left and 
right side of these two equations, we can get three equa- 
tions corresponding to Eqs. (Al), (A2) and (A3) 



exp [2AL] 



' 4a 2 + (A — c) A — c 2A — c 
4a 2 + (A + c) 2 A + c2A + c' 



(A4) 



exp 



2i arctan 



2a 



= exp 

exp [i2aL] 
= exp 



/ A+c A - c 

i arctan — arctan — — 

V 2a 2a 



„ / A + c A-c 

2i arctan — arctan — — 

V 2a 2a 



(A5) 
.(A6) 



Observing that there are two unknown parameters a and 
A but they need fulfill three equations, there are no so- 
lutions to the above equations for an arbitrary c. This 
means that the three-string assumption is not a solution 
to the Bethe ansatz equations. Also, we have solved the 
above equations numerically, but no solutions are found. 

Next we assume that the solution for a 3-atom system 
has the following form k\ = a — iA, ki = a + iA and 
k 3 = [3 with A > 0. With similar procedure as above, 
the Bethe ansatz equations can be represented as 



(2a + (3)L = (2n + n r )n + 2 ( 
A-c 



A + c 

arctan - 

\ a + 13 



— arctan 



, arctan — ) , (A7) 
a + (3 2a J ' v ' 



(2a - 13) L = (2n - n') ir + 2 ^arctan 



A 



— arctan 



A-c 

a- (3 



a — (3 

-) 

2a J 



arctan — ) , (A8) 



exp[2A£] = 



(a + (3f + (A - cY 



1/2 



(a + pf + (A + cf 



(a - pf + (A- c)' 



1/2 



1/2 



2A 



f R\ 2 ^(K J_ ^2l 1 / 2 2A + c 
(a- [3) + (A + c) 



(A9) 



where n and n' are integers which can be determined in 
the limit c — > 0. Here (2a + [3) corresponds to the total 
momentum which is conserved in the periodic boundary 
condition case, but not in the open boundary case. We 
will see that n = n' = 1 corresponds to the ground state. 
For the limit c = 0, we have 2a + (3 = 3tt/L which means 
that three atoms all occupy the lowest energy level. In 
the other limit c — > — oo, it is not obvious. By numeri- 
cally solving the Bethe ansatz equations, we can get the 
values of a, (3 and A. Corresponding to n = n' = 1, we 
get three sets of solution as shown in Fig. 3. Among the 
three solutions, we get 2a + (3 — ir/L corresponding to 
the solution A = — c + S (the solution denoted by squares 
in Fig. 3), where S is a small number in the large c limit. 
This solution is the trimer-like solution describing the 
ground state of the attractive 3-atom system. 
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